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EFFICIENT APPROXIMATION OF THE SOLUTION OF 
CERTAIN NONLINEAR REACTION DIFFUSION EQUATION II: 
THE CASE OF LARGE ABSORPTION 

EZEQUIEL DRATMANi 



Abstract. We study the positive stationary solutions of a standard finite- 
difference discretization of the semiUnear heat equation with nonlinear Neu- 
mann boundary conditions. We prove that, if the absorption is large enough, 
^Ji. ■ compared with the flux in the boundary, there exists a unique solution of such 

' a discretization, which approximates the unique positive stationary solution of 

1 ^^ ^ ' the "continuous" equation. Furthermore, we exhibit an algorithm computing 

an e-approximation of such a solution by means of a homotopy continuation 
method. The cost of our algorithm is linear in the number of nodes involved in 
' the discretization and the logarithm of the number of digits of approximation 

required. 
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1. Introduction 

This article deals with the following semilinear heat equation with Neumann 
boundary conditions: 

Ut = u^x-g\{u) in (0, 1) X [0,r), 

u^{\,t) = ag^{u{\,t)) in[0,r), 
u^(0,t) = in[0,r), 
u(a;,0) = uo(x)>0 in [0,1], 



where 51,52 S C'^(]R) are analytic functions in 2; = and a is a positive constant. 
The nonlinear heat equation models many physical, biological and engineering phe- 
nomena, such as heat conduction (see, e.g., |Can84[ §20.3], |Pao921 §1.1]), chem- 
ical reactions and combustion (see, e.g., [BE89| §5.5], jGri96[ §1-7]), growth and 
. migration of populations (see, e.g., [Mur02, Chapter 13], jPao92| §1.1]), etc. In 

particular, "power-law" nonlinearities have long been of interest as a tractable pro- 
totype of general polynomial nonlinearities (see, e.g., |BE891 §5.5], |GK041 Chapter 
7], |Lev90) . |SGKM95j . |Pao92l §1.1]). 

The long-time behavior of the solutions of ((T|) has been intensively studied (see, 
e.g., ICFQ91I , |GMW93) . |Qui93| , }R,os98) . [FROl] . [KTOl] . |AMTR,02] . [CQOi] and 
the references therein). In order to describe the dynamic behavior of the solutions 
of (HI) it is usually necessary to analyze the behavior of the corresponding stationary 



1991 Mathematics Subject Classification. 65H10, 65L10, 65L12, 65H20, 65Y20. 

Key words and phrases. Two-point boundary-value problem, finite differences, Neumann 
boundary condition, stationary solution, homotopy continuation, polynomial system solving, con- 
dition number, complexity. 

Research was partially supported by the following grants; UNGS 30/3084, UNGS 30/1066, 
CIC (2007-2009) and PIP 11220090100421. 

1 



2 



E. DRATMAN 



solutions (see, e.g., |FR01| . |CFQ91| ), i.e., the positive solutions of the following 
two-point boundary-value problem: 

{Uxx = 9i{u) in (0,1), 

^^.(1) - «52(w(l)), 
u,(0) = 0. 

The usual numerical approach to the solution of ([T]) consists of considering a 
second-order finite-difference discretization in the variable x, with a uniform mesh, 
keeping the variable t continuous (see, e.g., |BB98) ). This semi-discretization in 
space leads to the following initial- value problem: 

■^{U2 - Ui) ~ .gi(Ml), 

■^{uk+i - 2uk + Uk-i) " gi(wfc), (2 < fc < n - 1) 

^{Un-l - U„) - .gi(Wn) + ^92{u„), 

uoixk), {1 <k <n) 

where h :— l/{n — 1) and xi,...,Xn define a uniform partition of the interval 
[0, 1]. A similar analysis to that in |DM09| shows the convergence of the positive 
solutions of ^ to those of ([T]) and proves that every bounded solution of ^ tends 
to a stationary solution of ([3]), namely to a solution of 

^(u2 - Ml) - .gi(wi), 

^(ufc+1 - 2uk + Uk-i) - giiuk), (2 < fc < n - 1) 

■^{Un^l - Un) - gi{Un) + ^g2(Wn)- 

Hence, the dynamic behavior of the positive solutions of ([3]) is rather determined 
by the set of solutions (ui, . . . , u„) S (IR->o)" of 

Very little is known concerning the study of the stationary solutions of ^ and the 
comparison between the stationary solutions of ([3]) and ([1|). In IFROl] . [DM09j and 
[DralO] there is a complete study of the positive solutions of @ for the particular 
case gi{x) :— and 52(2:) :— x^, i.e., a complete study of the positive solutions of 

{0 = ■p-(M2 - ui) - w?, 
= ^{uk+i-2uk + Uk-i)-ul, (2<fc<n-l) 
= ^(u„_i-7.„)-< + f <. 

In |FR01j it is shown that there are spurious solutions of (|4]) for g < p < 2q — 1, 
that is, positive solutions of ^ not converging to any solution of © as the mesh 
size h tends to zero. 

In |DM09| and |DralO) there is a complete study of (|4]) for p > 2g — 1 and p < q. 
In these articles it is shown that in such cases there exists exactly one positive 
real solution. Furthermore, a numeric algorithm solving a given instance of the 
problem under consideration with n'^^^^ operations is proposed. In particular, the 
algorithm of [DralOj has linear cost in n, that is, this algorithm gives a numerical 
approximation of the desired solution with 0{n) operations. 

We observe that the family of systems ([5]) has typically an exponential num- 
ber 0(p") of complex solutions ( |DDM05j ). and hence it is ill conditioned from the 
point of view of its solution by the so-called robust universal algorithms (cf. [ParOO] , 
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[CGH+03) . [DMW09| ). An example of such algorithms is that of general continu- 
ation methods (see, e.g., jAG90| ). This shows the need of algorithms specifically 
designed to compute positive solutions of "structured" systems like 

Continuation methods aimed at approximating the real solutions of nonlinear 
systems arising from a discretization of two-point boundary-value problems for 
second-order ordinary differential equations have been considered in the literature 
(see, e.g., jABSWOB] . |Duv90| . jWatSOj ). These works are usually concerned with 
Dirichlet problems involving an equation of the form u^x = f{x, u, Ux) for which the 
existence and uniqueness of solutions is known. Further, they focus on the existence 
of a suitable homotopy path rather on the cost of the underlying continuation 
algorithm. As a consequence, they do not seem to be suitable for the solution 
of 11]). On the other hand, it is worth mentioning the analysis of [Kac02| on the 
complexity of shooting methods for two-point boundary value problems. 

Let (71,52 G C"^(R) be analytic functions in x = such that gi{Q) = 0, g[{x) > 0, 
g'/{x) > and g'/'ix) > for all x > with i = 1,2. We observe that gi and 
52 are a wide generalization of the monomial functions of system ([5]). Moreover, 
we shall assume throughout the paper that the functions g := (71/52 and G := 
Gi/gl are strictly increasing, where Gi is the primitive function of gi such that 
Gi(0) = 0, generalizing thus the relation 2q ~ 1 < p in In this article we 
study the existence and uniqueness of the positive solutions of (jH), and we obtain 
numerical approximations of these solutions using homotopy methods. In [?] there 
is a complete study of (jj]) for g := (71/(72 strictly decreasing. Furthermore, a similar 
analysis to that in [FROlj shows that there are spurious solutions of ^ for g 
strictly increasing and G strictly decreasing; i.e., the generalization of the relations 
g<p<2(7— lin®. According to these remarks, we have a complete outlook 
about the existence and uniqueness of the positive solutions of Q. 

1.1. Our contributions. In the first part of the article we prove that ^ has a 
unique positive solution, and we obtain upper and lower bounds for this solution 
independents of h, generahzing the results of [DM09) . 

In the second part of the article we exhibit an algorithm which computes an 
£-approximation of the positive solution of Q . Such an algorithm is a continuation 
method that tracks the positive real path determined by the smooth homotopy 
obtained by considering (jH) as a family of systems parametrized by a. Its cost is 
roughly of n log log e arithmetic operations, improving thus the exponential cost of 
general continuation methods. 

The cost estimate of our algorithm is based on an analysis of the condition num- 
ber of the corresponding homotopy path, which might be of independent interest. 
We prove that such a condition number can be bounded by a quantity indepen- 
dent of h := 1/n. This in particular implies that each member of the family of 
systems under consideration is significantly better conditioned than both an "aver- 
age" dense system (see, e.g., |BCSS98| Chapter 13, Theorem 1]) and an "average" 
sparse system ( |MR04[ Theorem 1]). 

1.2. Outline of the paper. In Section [2?2] we obtain upper and lower bounds for 
the coordinates of the positive solution of (jlj. 

Section [2731 is devoted to determine the number of positive solutions of (jlj. For 
this purpose, we prove that the homotopy of systems mentioned above is smooth 
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(Theorem [TH) . From this result we deduce the existence and uniqueness of the 
positive solutions of (H]). 

In Section 13] we obtain estimates on the condition number of the homotopy path 
considered in the previous section (Theorem I20[) . Such estimates are applied in 
Section |4] in order to estimate the cost of the homotopy continuation method for 
computing the positive solution of (j4|. 

2. Existence and uniqueness of stationary solutions 

Let Ul, . . . ,Un he indeterminates over R. Let gi and 52 be two functions of class 
C^(R) such that .gi(0) = 0, g'^ix) > 0, (x) > y g'/'ix) > for ah a; > with 
1 = 1,2. As stated in the introduction, we are interested in the positive solutions of 
(HI for a given positive value of a, that is, in the positive solutions of the nonlinear 
system 

(6) I = -{U^+i-2Uk + Uk-i) + h''gi{Uk), (2 < fc < n - 1) 

[ = -{Un-i-Un)!^gi{Un)-hag2{Un), 

for a given value a — a* > 0, where h :— l/{n — 1). Observe that, as a runs 
through all possible values in M>o, one may consider ([5]) as a family of nonlinear 
systems parametrized by a, namely, 

[ = -{U2~Ui) + ^giiUi), 

(7) < = -{Uk+i-2Uk + Uk-i) + h^gi{Uk), (2 < fc < n - 1) 

[ = -{Un-i-U„)!^gi{Un)-hAg2{Un), 
where A is a new indeterminate. 

2.1. Preliminary analysis. Let A, Ui, . . . , t/„ be indeterminates over R, set U :— 
{Ul, . . . ,Un) and denote by F : W^-^ — )• R" the nonlinear map defined by the 
right-hand side of ([T]). From the first n — 1 equations of (O we easily see that, for 
a given positive value Ui — ui, the (positive) values of U2, ■ ■ ■ ,Un, A are uniquely 
determined. Therefore, letting Ui vary, we may consider U2, ■ ■ ■ , Un, A as functions 
of Ul, which are indeed recursively defined as follows: 

Ui{ui) := Ul, 

t/2(wi) := ui + ^gi{ui), 

Uk+iiui) := 2Ukiui)-Uk-iiui) + h^gi{Ukiui)), (2<fc<n-l), 

Aim) := (^j^iUn-Un-iKui) + ^gi{Uniui))^/g2{Uniui)). 
Arguing recursively, one deduces the following lemma (cf. |DDM05l Remark 20]). 
Lemma 1. For any ui > 0, the following assertions hold: 

(t) [Uk - Uk-i){ui) - h^{yiiui) + J:'Z2 9i{Ujiui))) > 0, 
(it) Uk{ui)=Ui + h^(^^gi{ui) + Y!]ll{k-j)gi{Uj{ui))) > 0, 
(m) {Ul - Ul_,){ui) = hWg[{ui) + Y^^l 9'i{U,{ui))U'^{ui)) > 0, 
{^v) Ul{ui) = 1 + h\!^g'i{ui) + Y^U^k - j)g'i{U,{ui))U'^{ui)) > 1, 
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/or 2 < fc < n. 

As in [?] we have the fohowing lemma. This result is important for the existence 
and uniqueness of the positive solutions of ([T]). 

Lemma 2. For any ui > 0, the following assertions hold: 

(-) (t(i&)V)<o, 
M (%^)'(«i)>o, 

(-) (f7WrT)'K)>0' 

/or 2 < k < n. 

The next result studies the monotony of certain relations between gi and g2- 

Lemma 3. Let Gi he the primitive function of gi such that 6*1(0) = 0. If x > 0, 
then: 

W (^)V)>o. 

(m) //(f)'(x)>0, then (|i)'(a;)>0. 

(Hi) If exists d e [0,1) such that { \n{Gi{x) / gl{x)))' > 0, then (^^^ (x) > 0. 

Proof. Since gi is a positive and strictly convex function in ]R>o and .91 (0) = 0, we 
have that 

(9) ^ = r g,{t)g[{t)dt < g[{x) f g,{i)dt = g'i(x)Gi(x). 

Multiplying both sides by gi{x), we obtain 

giW 2gi(x)g^(a;) 
Gi(a;) 5?(x) ' 

which proves (0. 

Now suppose that {Gi/ gD' [x) > 0, then 

252(a;)52(2;) ^ 5i(a;) 



Combining this inequality with we obtain 

92ix) ^ 9l{x) g'ljx) ^ g'ljx) 
g2{x) 2Gi{x)g[{x) gi{x) ~ gi{x)' 

and is proved. 

Finally, suppose that exists d G [0,1) such that (ln(Gf (a:;)/g|(a;))) > 0, we 
deduce 

0<{HGiix)/gUx))y = d{HG^{x))y - (HgUx)))' 

^ ' ^ (ln(Gi(x)))'^ 
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Since (ln(Gi(x)))' — gi{x)/Gi{x) > 0, we have that 



(lii(GiW)) 

From this inequality, we obtain 

(10) '-^^^ ^ {H9im' < (ln(G,(.)))' = Me) , 

which completes the proof. □ 

2.1.1. Analogy between discrete and continuous solutions. Set Uk ■= Uk{ui) for 
2 < k < n. The first step in our analysis of the positive solutions of ([7]) is to 
estimate the discrete derivative of the solution u of Multiplying the identity 
u" = gi(u) by u' and integrating over the interval [0, x] it follows that 

(11) -u'{xf= u'{s)u"{s)ds^ u'{s)gi{u{s))ds^Gi{u{x))-Gi{u{0)) 
2 Jo Jo 

holds for any x G (0, 1), where Gi is the primitive function of gi such that Gi(0) — 0. 

The following result shows that ^( """^""^ )^, a discretization of ^u'{x)'^, equals 

the trapezoidal rule applied to gi{u{s))u' {s)ds up to a certain error term. 

Lemma 4. For every ui > and every 2 < m < n, we have: 

lfUm-Um-i\'^ gi{uk+i) + giiuk) , , 

(12) 2(^^) - g 2 

gi(m) r X 5i(um)/ N 

(U2 - Ui) (Um - Um-lj- 

Proof. Fix wi > and 2 < m < n. For m — 2 the statement holds by ([8]). Next 
suppose that m> 2 holds. From ([8|), we deduce the following identities 

2(m2 -ui) 

^^2 9i{ui), 

^ 9i(Uk), (2<fc<TO-l). 

Multiplying the first identity by {u2 — ui)/Ah and the fcth identity by {uk+i — 
Uk^i)/2h, we obtain 

1 / /Uk+l-Uk\^ /Uk~Uk-l\^\ IfUk+l-Uk Uk-Uk-l\ , . 

for 2 < k < m ^ \. Note that {uk+i — Uk-i)/2h is a numerical approximation of 
u'{{k— l)h) for 2 < fc < TO — 1. Adding these identities multiplied by h we obtain: 

l/u™-'u„_i\2 gi{ui) ^ , V^\9l("fc)/ , ^ 

^(^ ) = {U2-Ul)+ 2^ ^ [Uk+l-Uk+Uk- Uk-l) 

k=2 

'v" .9i("fc) +.9i(ufc+i) . . 9i{ui) , s 

= 2^ ^ (Wfc+l-Ufe) —{U2-Ui) 

fc=l 

^ ("m-Mm-l)- 

This finishes the proof of the lemma. □ 
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Substituting x for 1 in wc obtain the foUowing identity (cf. |CFQ91[ §3]): 
(13) \a^gl{u{l)) = G,{u{l))-G,{um. 

From this identity one easily deduces that u(l) is determined in terms of I'o '■— u(0), 
say, u(l) = /(t'o)- Therefore, by it is possible to restate ^ as an initial-value 
problem, namely 

Uxx = 9i{u) in (0, 1), 
u(0) ^ vq, 
u,(0) = 0, 

where i^o > is the solution of the equation Ux{^) = ag2{f{vQ))- Our purpose is to 
obtain a discrete analogue of this identity, which will be crucial to determine the 
values Ml of the positive solutions of (O . 

Let (a, u) := (a, ui, . . . , u„) G (R>o)"^^ be a solution of ([7]). From the last 
equation of ([7]) we obtain 

l/u„-w„_i\2 If , s h , .\'^ 



/U„-W„_l\2 1/ , X I \\' 

\ = 2(«.92(-«)-25iK)j 

2(U„) + -;^g\{Un){^-g\{Un)-ag2{Un)^-—g\{Ur,) 
I99/ \ ^ / \ ( ^n — 1 \ ^ 9 / \ 

Combining this identity with Lemma |4] we obtain 

1 2 2/ N 51K) +.9l("fc+l)/ ^ .91 (ui). h'^ 2, N 

-a ffaK) = 2^ ^ (Mfe+i - Uk) - ui) + — 5i(u„). 

fc=i 

Using the identity gi{ui)[u2 — ui) = ^^^(mi), we deduce that 

1 h"^ 

-a^gliun) - (Gi(u„) - Gi(ui)) = S + —{gl{un) - 5i(ui)), 

holds, where Gi is the primitive function of gi such that Gi(0) = and E is defined 
as follow: 

(14) ^ 2^ 2 \'^k+'^ - ""fc) - (Gi(u„) - Gi(Mi)j. 

It is easy to check that E is the error of the approximation by the trapezoidal rule 
of the integral of the function gi in the interval [ui, w„], considering the subdivision 
of [ui, Un] defined by the nodes wi, . . . , u„. Moreover, taking into account that gi is 
a convex function in ]R>o, we easily conclude that E >Q holds. Therefore, from the 
previous considerations we deduce the following proposition, which is the discrete 
version of (fT3l). 



Proposition 5. Let [a.u) E (K>o)"+^ be a solution of Then 

(15) la^gUun) - (Gi(m„) - Gi(ui)) =E+^{gl{ur.) - gUui)), 

where Gi is the primitive function of gi such that Gi(0) = 0, and E is defined as 
in (|14p . Furthermore, if we consider E as a function of ui according to (jl4p . where 
Uk ■— Uk{ui) is defined as in ([8]) for 2 < k < n, then E is a positive increasing 
function over R>0. 
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Proof. For the above considerations, we only have to prove that E is an increasing 
function over M>0. 

We consider i? as a function of ui, where Uk Uk(ui) is defined as in (jS]) for 
2 < k < n. If rewrite E as follows 

(16) E = g _ _ 

fe=i 

then it suffices to show that each term of the previous sum is an increasing function 
over R>0. In fact, fixl<fc<n — 1; the derivative of the fcth term of (|T6| as 
function of ui is 

d f gi{uk) + gi{uk+i) , \ (r< f \ n t \\\ 
^ [uk+i - Uk) - (Gi(ufc+i) -Gi[uk))j = 

9i{uk)v'k + g'i{uk+i)v'k+-^ .gi(ufe) +.gi(ufe+i) . , ,^ 
= ^ K+i - ""fc) + ^ K+i - 

-(gi(wfe+iX+i - gi(wfe)w;), 

where := U'f.{ui) and f^^^ := C/^^]^(ui). Adding and subtracting v'f. in each 
occurrence of w^^-^ we obtain 

d [ gi{uk) + giiuk+i) . \ (n t \ n i \\\ 
^ (ufc+i -ufe) - (Gi(ufe+i) - Gi(ufe)) j = 

( g'liuk) ^ g'i(uk+i) , \ ( I \ I \\\ I 

= y ^ (wfe+i - ufc) - (5i(ufc+i) - gi{uk)) jWfc 

, /ffUwfe+i)/ N gi{uk+i) - gi{uk)\, , 

+ (^ ^ {Uk+l-Uk) ){Vk+l-Vk) 

( g'Mk) ^ g'i{uk+i) , \ ( / \ I \\\ I 

= \ ^ - Uk) - [gi(uk+i) - giiuk)) jwfc 

, ( g'liuk+i) g'i{£,k) \i w , ,^ 

+ (^ 2 —){Uk+l -UA;)(Wfc+l -^'fc), 

where ^k G (wfc,Wfe+i) is obtained after applying the Mean Value Theorem to 
gi{uk+i) — gi{uk)- It is easy to check that 

g'ljuk) + g'ljuk+i) f . ( , . / ^^ 
[Uk+i - Uk) - [giiuk+i) - gi(Uk)) 

is the error of the approximation by the trapezoidal rule of the integral of the 
function g[ in the interval [uk, Uk+i], and the convexity of g'l ensures their positivity. 
In the other hand, since g'l is increasing, we have that gi{uk+i) — gi{S,k) > 0. Finally, 
from Lema[Tl {v'i^_^-^ — v'f.), (uk+i — Uk) and v'j^ are positive numbers for ui € M>0. 
Therefore, the fcth term of p6p is increasing over ]R>0 for 1 < fc < n — 1, which 
completes the proof. □ 

2.2. Bounds for the positive solutions. In this section we show bounds for the 
positive solutions of ^ . More precisely, we find an interval containing the positive 
solutions of ([7]) whose endpoints only depend on a. These bounds will allow us to 
establish an efficient procedure of approximation of this solution. 
Let g : IR>o — IR.>o and G : M>o IR>o be the functions defined by 

(17) g{x) := ^(x), 

52 



9 



and 

(18) G(x) 

where Gi is the primitive function of gi such that Gi(0) = 0. 
As in [?, Lemma 7] we have the fonowing resuh 

Lemma 6. Let (a, u) E (M>o)"^^ be a solution of ([7]) for A = a. Then 

ag2{u„) < gi[un). 

From Lemma IHl we obtain the following corollary. 

Corollary 7. Let (a, u) e (M>o)"+^ he a solution of ([7]) for A — a. If the function 
g defined in \n\j is surjective and strictly increasing, then 

Let (a, u) e (R>o)"^"'^ be a solution of ^ ioT A — a. As in [?, Lemma 9] we 
obtain an upper bound of Un in terms of ui and a. 

Lemma 8. Let {a, u) £ (R>o)"^^ be a solution of ([7]) for A — a, and let C{a) be 
an upper bound o/u„. Then u„ < e^^ui holds, with M :— g[(^C{a)^. 

The next lemma shows a lower bound of ui in terms of a. 

Lemma 9. Let {a, u) e (R>o)"^^ be a solution of ^ for A — a, and let C{a) be 
an upper bound of Un ■ If the function g defined in |_?7p is surjective and strictly 
increasing, then 

holds, where M := g[(^C{a)^. 

Proof. From Lemma |8] and Corollary [7] we deduce 

g^^ia) < Un < e^'ui, 
which immediately implies the statement of the lema. □ 

In the next lemma we obtain another upper bound of u„ in terms of ui and a. 
This upper bound will allow us to find an upper bound of u„ in terms of a. 

Lemma 10. Let {a, u) e (R>o)"^^ be a solution of ([7]) for A — a. Then 

G(u„) < G(ui) + —, 

where G is defined in hl8\). 

Moreover, if G is surjective and strictly increasing, then 

2 

u„ <g-i(gk) + ^). 

Proof. From Proposicin[3]and Lema[T] we deduce the following inequality 

Gi(u„) - —gl{un) < Gi(ui), 

where Gi is the primitive function of gi such that Gi(0) = 0. Dividing for g2(un), 
we obtain 

Gi(iti) 
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Since (?| is an increasing function, we conclude that 



2 gl{un) 



which prove the first part of the lemma. 

Now, suppose that G is surjective and strictly increasing, then G is an invertible 
function and their inverse is strictly increasing. Combining this remark with the 
last inequality, we obtain 

2 

u„<G-i(GK) + y), 
and the proof is complete. □ 

From this lemma we obtain upper bounds for u\ and u„ in terms of a. 

Proposition 11. Let (q!,u) G (M>o)"+^ he a solution of ^ and let g and G be the 
functions defined in ^17^ and il8\) respectively. Suppose that 

• exists d e [0, 1) such that (ln{Gf{x) / gjix)))' > for all x > 0, 

• G"{x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then 



1 - d 

Moreover, if g and G are surjective functions, then 



Proof. Combining Lemma |H] and the Mean Value Theorem, we deduce that exists 
^ > between ui and u„ such that 

G'(OK - "1) = GK) - G(ui) < 
By Lemma [1T|m|1 . we have that 

1 n—l I N 

(Un-ui) = h (^—^31(^1) + 2^[n- j)gi{uj) j > — — > 0. 
Combining both inequalities, we obtain that 



G'(0^<G'(OK-.i)<^. 



11 



Since G"{x) > for all a; > 0, we see that G' is an increasing function. Furthermore, 
we have the following inequality: 

2^^,^ N i ^ / gi(m)g|("i) ~ gl(^l)2g2(^l)g2(^l) ^ f ^ 

V 52(^1) ^ 

/ G'i(Mi)2g2(Mi)g2(Mi) \ .9i(m) 
^^^> V .giK).9iK) ^9i(«i) 

/ (ln(.9|K)))\g^(ui) 



(21) a^>{l-d)^i^ = {l-d)gHu,), 



(ln(Gi(wi)))'^5i("i)' 
Taking into account the first condition of the statement, we deduce that 
0<{HGiix)/gUx))y = d(ln(Gi(x)))'-(ln(gi(x)))' 

= (ln(G,(.)))'(rf-lM^ 

Since (ln(Gi(a;))) = gi(a;)/Gi(x) > 0, we conclude that 

20 ^ ^^^^ , < d. 

(HGiix)))' 

Combining ([T9l) and ([20]) . we obtain 

j .9i(m) 
gi(Mi) 

and the first assertion of the proposition is proved. 

Now, suppose that g and G are surjective functions. From ([3]), we have that 
g and G are strictly increasing functions. Combining this remark with (|2ip and 
Lemma [TOl we obtain the desired upper bounds for ui and u„. □ 

Combining Proposition 1111 and Lemma [9] we obtain the following result. 

Lemma 12. Let (a, u) G (R>o)"^^ &e a solution of ([7]) and /ei g and G be the 
functions defined in ^17^ and il8\) respectively. Suppose that 

• G and g are surjective functions, 

• exists d e [0, 1) such that ( ln(Gf > for all x > 0, 

• G"{x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then 

ui < g~'^{aC{a)), 

where 

_ g2(Ci(a))a2 
^^"^ ' ^ ^ 2g2(5-K«)/e^0G'(g-K")/e^0' 

an^^ A'f g[{Ci{a)). Furthermore, 

2 

<G-i(G(.g-i(aG(a))) +^). 
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Proof. Let {a,u) G (R>o)"^^ be a solution of ([7]). From Lemmas IT] and ITOl we 
deduce the inequalities 

• 5i(ui) < ^(^5i("i) + .9i("2) H ^ gi{un-i) + ^gi{un)j = ag2(u«), 

• un < G-i(G(iii) + 4)- 
Combining both inequalities we obtain 

, , ff2K) 32(g-(g(.o + 4 

g(ui) < a — — r- < a 



52(wi) ff2(til) 

Since (72 is an increasing function in M>o and (G^^)' is a decreasing function in 
by the Mean Value Theorem, we obtain the following estimates 

32(^1) + .g^(G-i(G(ui) + ^^)) (g-i(gK) + 4) - "1 



g{ui) < a 



< a 



< a 



52(^1) 

32^) + g^(G-i(G(ui) + ^))(G-i)'(GK))- 



52 ("1) 

g^(G-i(GK) + 4))^ 



2g2{u^)G'{ui) 
From Proposition [TT] and Lemma [9] we conclude that 

Ml < 5"^(aG(a)), 

where 



G{a) := 1 



2g2(5~i(a)/e^0G'(.g~i(«)/e*0' 
with M := 5i(Gi(a)) and 



Gi(a) :=G-i(G(g-i(- 



Combining this remark with Lemma [TU] we obtain 



Un < G 



''[G[g''[ad{a))) 



a" 



2 

which immediately implies the statement of the lemma. □ 

2.3. Existence and uniqueness. Let P : (K>o)^ — >■ R be the nonlinear map 
defined by 

(22) P(a,ui) i(C/„_i(ui) - Un{ui)) ~ §51 (C^nK)) + "52 ([/^(mi)) • 

Observe that U\) — ^ represents the minimal equation satisfied by the coor- 
dinates (a, Ul) of any (complex) solution of the nonlinear system ([7]). Therefore, 
for fixed a G M>o, the positive roots of Pia, U\) are the values of u\ we want to 
obtain. Furthermore, from the parametrizations ([5]) of the coordinates 1*2, . . . ,u„ 
of a given solution (q;,ui, . . . ,m„) € (R>o)""^^ of ([7]) in terms of ui, we conclude 
that the number of positive roots of P(q;,C/i) determines the number of positive 
solutions of (O for such a value of a. 

Since P(^, C/i) is a continuous function in (R>o)^, as in [?, Proposition 4] we 
have the following result: 
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Proposition 13. Fix a > and n e N. // the function g defined in JT?! ) is 
surjective, then ^ has a positive solution with A = a. 

In order to establish the uniqueness, we prove that the homotopy path that we 
obtain by moving the parameter a in R>o is smooth. For this purpose, we show 
that the rational function A{Ui) implicitly defined by the equation P{A,Ui) — is 
increasing. We observe that an explicit expression for this function in terms oi Ui 
is obtained in ([5]). 

Theorem 14. Let A > Q he a given constant and let A[Ui) he the rational function 
of l^i). Let g and G he the functions defined in {11^ and i f_?($|) respectively. Suppose 
that 

• G and g are surjective functions, 

• exists d e [0, 1) such that {\n{G'l{x) / gl{x)))' > for all x > 0, 

• G"{x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there exists 
M{A) > such that the condition A' {ui) > is satisfied for n > 1 + M{A)/{2-2d) 
and ui e A-^ ((0, A^ n R>o. 

Proof. Let Ui,U2, ■ . ■ ,Un, A be the functions defined in ([S]). For ui > 0, we 
denote by /(wi) := G'i(L/„(ui)) — Gi(C/i(mi)) the integral of the function gi in 
[Ui{ui), Uniui)], and by T{ui) the trapezoidal rule applied to /(ui), with the nodes 
Ui{ui), t/2 (ui ),,... , Un{ui). More precisely, T is define as follows: 

rr. 9l{Uk + l)+gi{Uk) f,, . 

J •= 2^ 2 [Uk+i - Uk). 

k=l 

Finally, set E := T — L. Combining Proposition [5] and the convexity of gi, we 
deduce that E > and E' > in K>o, where E' represent the derivative of E with 
respect of ui. 

According to Proposition [Sj Ui,U2, ■ . ■ ,Un, A satisfy the discrete version ((T5|) of 
the energy conservation law pT|) . Dividing both sides of (|T5|) by G'i(C/„) we obtain 
the following identities: 

(23) ^^^2 9UUn) _ T ^ h^ gfjU^) / gf{U^ 



2" Gl{Un) Gl{Un) ' 8Gl(t/„)V^ gl{Un) 

Taking derivatives with respect to C/i at both sides of (P^. we have 



AA' 



2 VGi( 



.24^ Gi(C/„) 2 VGi([/„) 

^ ' _( T \' ( 9l{Un) \'( gl{Ui) \ gWn) ( gliUi) 
\GiiUn)) 8 \Gi{Un)) \ gliUn)) 8 Gi{U^)\gl{U^), 

Let ui e A"^((0,yl]) nR>o. By Lema[Tl C/j(ui) and U[{ui) are positive for 1 < i < 
n. Furthermore, 51, (72, g, G and Gi are positive and increasing functions in R>o. 
Throughout the demonstration we will use these conditions repeatedly. 

From Lemmas |3] and [2] we deduce that gJ(C/„)/Gi([/„) is an increasing function 
and gf{Ui)/gf{Un) is a decreasing function. Combining these remarks with ([M)) 
we obtain 

(-ig)'-»((s;^)'-T(lti)>""- 
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we see that the incquaUty above may be rewritten in the form 

(-asx- > (^(f)'-iii(^)>-' 

-(T(lti)')<-'- 

We claim that {T / Gi{Ui)) [ui) > 0. Indeed, by Lemmas [T] and [H we have that 

Combining this result with Lemma [31 we conclude that 
Combining the claim above with p5|) we deduce that 

(-ig)<".'>(5io(T)'4(iti)>-)- 

In order to prove the positivity of A'(ui), we rewrite the right side of ((26)l . 

Ig2([/„)V/J 2 VGi(c/„)J J^^'^- 



V GjiUn) 2 G?(C/„) 
/ E'l 



^JI-EV A^gl{U^){G,m)' ( {9l{Un))'G,{Un) \\ 
G?(C/„) + 2G?(C/„) I 5l(^n)(Gi(t/n))Vy 

Since there exists d £ [0, 1) such that ( ln(Gf (a:)/(7|(x)))' > for all a; > 0, we have 
that 



\ (gi(t/n))'Gi(L/„) 



<?2([/„)(Gi(!7„)) 



7 > 1 - d. 



Furthermore, by the positivity of E'iui) and the definition of /(ui), we conclude 
that (^E' I — EI'){ui) > —E{ui){G{Un)) [uij- From these inequalities, we deduce 
that 

(Ml) > 



VG2(C/„)V/y 2 VGi({/„) 



2 VGi({/„)^ ^ 



> 

Combining these remarks with (|26p . we obtain 
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From ([T5|) . it follows that 



If we prove that the second term of the right side of (P?]) is positive, we obtain 
that 

(■-||g)<«.»(^^^-^)<«.»». 

which immediately implies the statement of the theorem. Therefore, it suffices to 
show that 

Since gi is an increasing function, we only need to show that 
1 d 1 d 



k=l 

where 



gi(t/fc+i) +.gi(t/fc) ^^, ^ 
J-k ■= 2 (Uk+i-Uk), 

Ek Tk — Ik, 

with Ik Gi{Uk+i) ~ Gi{Uk)- Note that, for ui > 0, /fc(wi) is the integral of gi 
in [[/fe(wi), J7fe-|_i(ui)], Tfe(ui) is the trapezoidal rule applied to Ik{ui) and Ek{ui) 
is the error of such approximation. 

In order to prove that (1 — d)T(ui)/2 — E{ui) > 0, we show that 

(29) ^-^Tk{ui) - Ekiui) > 

holds for 1 < fc < n — 1 . By [?] , we have that 

Ek(ui) < \^ [Uk+i-Uk) j[ui). 

From Lemma [1] and the monotonicity of g[ , we deduce that 

Ekim) < (^i^(L/fc+i-C/,)^)(«i) 

< (/.gif^ ^^(^-^-),+ -^-^^-) (^..i-^..))K) 
g'liUn 



= {h'-^Tk)M 
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Thus, we see that ([29|) is satisfied if the inequality 
(30) h < — ^ 

holds. From Lemma 1121 and the monotonicity of g'l, we deduce that there exists 
a constant M{A) > independent of h such that g'i{Un){ui) < M{A) for m G 
A~^((0,^]) nM>o- This shows that a sufficient condition for the fulfillment of (1501) . 
and thus of A'{Ui) > 0, is that n-1 > M{A)/(2 - 2d) holds. This finishes the 
proof of the theorem. □ 

In order to prove the uniqueness of positive solutions of (O , we still need a result 
on the structure of the inverse image of A on the interval under consideration. 

Lemma 15. Let A> be a given constant and let A(Ui) be the rational function 
of (0^. Let g and G be the functions defined in ( [j7| j and respectively. Suppose 
that 

• G and g are surjective functions, 

• exists d e [0, 1) such that (ln(Gf (a:;)/gr|(a;)))' > for all x > 0, 

• G"{x) > for all x>0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there 
exists M{A) > that satisfies the following condition: for n > 1 + M(A)/{2 — 2d) 
there exi sts c c(n, A) > such that ^"^(0, A]) n M>o = (0, c] . 

Proof. By Theorem [T3] we have that there exists M (A) > such that the condition 
A'{ui) > is satisfied for n > 1 + MiA)/i2-2d) and ui e A-^ {{0, A]) CiR^o- Fix 
n > 1 + AI (A)/ {2 — 2d). From Lemma [1] we deduce that Un{ui) defines a bijective 
function in M>o and that 

^("i) = [l^9i{Ui{ui)) + J2 hgi{Uk{ui)) + ^giiUn{ui)))/g2iUM). 

k=2 

Since < Ui{ui) < ■ ■ ■ < C/„(ui), we have the following inequalities: 
^9{Un{ui)) < A{ui) < g{Uniui)). 

Since lim„^^o+ 9{Un{ui)) — 0, there exists e > such that (0,e] C A^^((0,^]) n 
R>o. We claim that 

(0,co] = A-i((o,^])nR>o 

with 

Co :=sup{e: (0,e] C A-i((0,^])nR>o}. 
Indeed, from the definition of cq we obtain that (0,co) C A~^((0,^]) n M>o and 
that lim^^_j^^- A{ui) — A{co) < A. Therefore, we deduce that 

(0,co] c A-i((o,^])nM>o. 

We now show that the last set inclusion is an equality. Suppose that there exists 
5 > Co such that A{S) < A. Let a inf{(5 : S > co,A{5) < A}. From the 
definition of cq, the interval (co,ci) is not empty. Since A{x) > A for all x G 
(co,ci), we have that A'{ci) < 0, which contradicts the fact that A'{ui) > for all 

wi e A-i((o,yt]) nR>o. □ 

Now we state and prove the main result of this section: 
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Theorem 16. Let a > be a given constant. Let g and G be the functions defined 
in |i7[ ) and I118\) respectively. Suppose that 

• G and g are surjective functions, 

• exists d e [0, 1) such that (lii(G'f (a;)/.g|(a;)))' > for all x > 0, 

• G"{x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there exists 
M(a) > such that Q has a unique positive solution for n > 1 + M{a)/{2 — 2d). 

Proof. Proposition [T51 shows that ^ has solutions in (M>o)" for any a > and 
any n S N. Therefore, there remains to show the uniqueness assertion. 

By Theorem [T^ we have that there exists M{A) > such that the condition 
A'{ui) > is satisfied for n > l + M{A)/{2-2d) and mi G ((0, y^]) nR>o. From 
Lemma [T51 there exists c — c{n, a) such that ^~^((0, a])nM>o = (0, c]. Arguing by 
contradiction, assume that there exist two distinct positive solutions (ui, . . . ,'u„), 
(ui, . . . , Un) & (IR>o)" of ([6]). This implies that ui ^ ui and A{ui) = A{ui), where 
A{Ui) is defined in ([8|). But this contradicts the fact that A'{ui) > holds in (0, c], 
showing thus the theorem. □ 



3. Numerical conditioning 

Let be given n S N and a* > 0. Let g and G be the functions defined in p7)) 
and (fTS|) respectively. Suppose that 

• G and g are surjective functions, 

• exists d e [0, 1) such that (ln(G'f (x)/g|(a:)))' > for all 2: > 0, 

• G"(x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi (0) = 0. In order 
to compute the positive solution of ([7]) for this value of n and A — a* , we shall 
consider ([7]) as a family of systems parametrized by the values a of A, following 
the positive real path determined by ([7]) when A runs through a suitable interval 
whose endpoints are a* and a*, where a* be a positive constant independent of h 
to be fixed in Section |4l 

A critical measure for the complexity of this procedure is the condition number of 
the path considered, which is essentially determined by the inverse of the Jacobian 
matrix of ^ with respect to the variables Ui, . . . ,Un, and the gradient vector of ([7]) 
with respect to the variable A on the path. In this section we prove the invertibility 
of such Jacobian matrix, and obtain an explicit form of its inverse. Then we obtain 
an upper bound on the condition number of the path under consideration. 

Let F :— F{A, U) : ]R"+^ — > M" be the nonlinear map defined by the right-hand 
side of ^ . In this section we analyze the invertibility of the Jacobian matrix of F 
with respect to the variables U , namely, 

ri -1 \ 
-1 '•. ■■. 
'■. ■•. -1 

with Fi := 1 + ^h'^g[{Ui), F, 2 + h'^g[{Ui) ioi 2 < i < n - I and F„ := 
l + ^h^g'iiUn)-hAg'^{U„). 
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We start relating the nonsingularity of the Jacobian matrix J{a, u) with that 
of the corresponding point in the path determined by Let {a,u) e (M>o)"^^ 
be a solution of ([7]) for A = a. Taking derivatives with respect to [/i in ([5]) and 
substituting ui for Ui we obtain the following tridiagonal system: 



-1 

V 



\ / 1 



-1 r„K) / 







\ hg2{U.n{ui))A'{ui) I 



For 1 < fc < n — 1, we denote by := Afc(A, U) the fcth principal minor of the 
matrix J [A, U), that is, the (fc x fc)-matrix formed by the first k rows and the first 
k columns of J {A, U). By the Cramer rule we deduce the identities: 



(31) hg2{Un{ui))A'{u,) 

(32) det(j(a,w))[/^(wi) 



det ( J(q!, u)) , 

hg2{Un{ui))A'{ui)det{Ak-i{a,u)), 



for 2<k<n. 

Let a > be a given constant. Then Theorem [T3] asserts that A'{ui) > 
holds. Combining this inequality with (I5T]) we conclude that det ( J(a, u)) > 
holds. Furthermore, by ([5^ . we have 



(33) 



C/fe(ui) = det (Afe_i(a,M)) (2 < fc < n). 



Combining Remark [T lfru)) and (l33l) it follows that det(Afc(a, u)) > holds for 
1 < k < n — 1. As a consequence, we have that all the principal minors of the 
symmetric matrix J(a, u) are positive. Then the Sylvester criterion shows that 
J(a,u) is positive definite. These remarks allows us to prove the following result. 

Theorem 17. Let {a,v) e (R>o)"^^ be a solution of ^ for A = a. Let g and G 
be the functions defined in ( [j7| ) and il8\) respectively. Suppose that 

• G and g are surjective functions, 

• exists d e [0, 1) such that ( ln(G'!^(a;)/.gf (a;)))' > for all x > 0, 

• G"{x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) — 0. Then there 
exists M{a) > such that the matrix J{a, u) is symmetric and positive definite for 
n > l + M{a)/{2 ~ 2d). 

Having shown the invertibility of the matrix J(a, u) for every solution (a, u) G 
(R>o)"'^^ of (O, the next step is to obtain explicitly the corresponding inverse 
matrices J^^{a,u). For this purpose, we establish a result on the structure of the 
matrix J^^{a, u). 

Proposition 18. Let (a, u) G (M>o)"^^ be a solution of Let g and G be the 

functions defined in Jj?] ) and il8\) respectively. Suppose that 

• G and g are surjective functions, 

• exists d G [0, 1) such that [\n{Gi{x) / gl{x)))' > for all x > 0, 

• G"{x) > for all x > 0, 
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hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there 
exists M(a) > such that the following matrix factorization holds: 



I !i2. 



1 



V 



1 



/ 



/ \ 1 "2 

/ \ d(J) d(.I) 



d{J) d{J) ■ ■ ■ d{J) d(J) / 

for n > 1 + M{a)/{2 - 2d), where d{J) := det (j{a,u)) and u'j. := U'f.{ui) for 
2<k<n. 

Proof. Since J{a,u) is symmetric, invertible, tridiagonal and their (n — l)th prin- 
cipal minor is positive definite, the proof follows in the same way as that of [DM09[ 
Proposition 25]. □ 

From the explicitation of the inverse of the Jacobian matrix J{A, U) on the 
points of the real path determined by ([7]), we can finally obtain estimates on the 
condition number of such a path. 

Let Of* > and a* > constants independents of h be given. Then Theorem 
1161 proves that ([7]) has a unique positive solution with A = a for every a in the 
real interval I := I(Q!*,a*) whose endpoints are a* and a*, which we denote by 
U2{ui{a)) , . . . , [/„(ui(q;))). We bound the condition number 

K := max{||v3'(Q;)||oo : a el}, 

associated to the function : I — > M", f{a) :— [ui{a), U2{ui{a)), . . . , [/„(ui(a))). 
For this purpose, from the Implicit Function Theorem we have 

dF 



J- 



We observe that {dF/dA){a,>f{a)) = (o, . . . ,0, -hg2{Un{ui{a))) 
Proposition [TH] we obtain 



holds. From 



hg2{Un{ui{a))) ( 

-T—rjj TT^ i ^' ^2 ("1 (a)) ,---,U„[Ui (a)) 

det [J[a, 'p(ce))) V ' ' 



ll'p'(a)lloo = 

Combining this identity with (|3ip . we conclude that 



ll^'(«)lloo = II A'iMa)) ^' ^^("^("))' ■ • ■ ' '^"(^1 
From Lemma [1] we deduce the following proposition. 

Proposition 19. Let a* > and a* > constants independents of h be given. Let 
g and G he the functions defined in ( [j7| j and H8\) respectively. Suppose that 

• G and g are surjective functions, 

• exists d e [0, 1) such that {\n{G'i{x) / gHx)))' > for all x > 0, 

• G"ix) > for all x > 0, 
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hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there 
exists M{T) > such that 

II ,1 ^|| U'n{ui{a)) 
A'[ui[a)) 

holds foraeIandn>l + M{I)/{2 - 2d). 
Combining Proposition [TH] and ([25)) we conclude that 

4Gi([/„(ui(a))) 



(/?'(q;)||oo < 



(1 - d)agi{Un{ui{a))) 
Applying Lemma |9] and Proposition [TT] we deduce the following result. 

Theorem 20. Let a* > and a* > constants independents of h he given. Let g 
and G be the functions defined in |_??| ) and kl^) respectively. Suppose that 

• G and g are surjective functions, 

• exists d € [0, 1) such that (ln(Gf > for all x > 0, 

• G"(x) > for all x > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then there 
exists a constant Ki(Q;*,a*) > independent of h such that 

K < ki(q!*, a*). 
4. An efficient numerical algorithm 

As a consequence of the well conditioning of the positive solutions of ([7]), we 
shall exhibit an algorithm computing the positive solution of ([7]) for A = a* . This 
algorithm is a homotopy continuation method (see, e.g., |OR70[ §10.4], |BCSS98l 
§14.3]) having a cost which is linear in n. 

There are two different approaches to estimate the cost of our procedure: using 
Kantorovich-type estimates as in |OR70[ §10.4], and using Smale-type estimates 
as in |BCSS98l §14.3]. We shall use the former, since we are able to control the 
condition number in suitable neighborhoods of the real paths determined by ([T]). 
Furthermore, the latter does not provide significantly better estimates. 

Let a, > be a constant independent of h. Let g and G be the functions defined 
in (fT7|) and (fT8|) respectively. Suppose that 

• G and g are surjective functions, 

• exists d G [0, 1) such that (ln(Gf (a;)/g|(a;)))' > for all a; > 0, 

• G"(x) > and g"{x) > for all a; > 0, 

hold, where Gi is the primitive function of gi such that Gi(0) = 0. Then the 
path defined by the positive solutions of (O with a € [a*, a*] is smooth, and the 
estimate of Theorem 1201 hold. Assume that we are given a suitable approximation 
of the positive solution (p{ar) of ([7]) for A = a,. In this section we exhibit an 
algorithm which, on input u^'^\ computes an approximation of ip{a*). We recall 
that ip denotes the function which maps each a > to the positive solution of ([7]) 
for A = a. More precisely ip : [a*, a*] — > R" is the function which maps each 
a e [a,,Q*] to the positive solution of ^ for A — a, namely 

(p{a) := (ui(a), . . . := {ui{a),U2{ui{a)) , . . . ,Un{ui{a))) . 

From Lemma [12] and Lemma [9l we have that the coordinates of the positive 
solution of ([7|) tend to zero when a tends to zero. Therefore, for a small enough. 



21 



we obtain a suitable approximation of the positive solution ^ for A — a^, and we 
track the positive real path determined by ([7]) until A — a* . 

Let < a* < a* be a constant independent of h to be determined. Fix a € 
[a»,a*]. By Lemma [121 it follows that ip{a) is an interior point of the compact set 



where 



with 



if„:={ueM":||u||oo<2C2(a)}, 



C{a) := 1 + 



9'2{Ci{a))c 



2g2(.9-H«)/e^0G'(.9-K")/e^0' 



a 
~2 



and Af := g[{Ciia)). 

First we prove that the Jacobian matrix Ja{u) :~ {dF / dU){a, u) is invertible in 
a suitable subset of ■ Let u G M" and v G R" be points with 

\\u- Lp{a)\\oo < 5a, ||u - '^(q;)||oo < (^Q, 

where 5q > is a constant to be determined. Note that if 5p < C2(a) then u G Ka 
and V S Ka ■ By the Mean Value Theorem, we see that the entries of the diagonal 
matrix Ja{u) — Ja{v) satisfy the estimates 

{Jaiu)- Ja{v))^\^ < 2//,V(2C72(«))<^a, (l<i<n-l) 

(Ja(u)-Ja («))„„ I < 2/imax{«5^'(2C2(a)),<(2a2(«))}<5„. 

By Theorem 1171 and Proposition [18] we have that the matrix J^p(a) ■— Jaiy^ict)) = 
{OF /dU){a, (p{a)) is invertible and 



E 



C// {ui (a)) f/j {ui (a)) U.l {m (a)) C/j (t/i (a)) 



{.J} det(J^(„)) 

holds for 1 < i,j < n. According to Lemma [TJ we have C/4(wi(a)) > 
U2{ui{a)^ > 1. These remarks show that 



> 



(34) 



< r]aSa 2 



< 



j:;=2h'U',iui{a))+hU:,{u^{a)) 
I det(Jy(Q))| 



< 2riaSa 1 



where ?7q 
identity: 



I det(Jy(„))|y' ' 

2max{<(2C2(a)),a9^'(2C2(a))}. From we obtain the following 



hU'^{u,{a)) _ 

|det(J<^(o))| A' {ui{a)) g2{un{a)) 
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From ([28l) . we have that 
(35) ''^UMc^)) _ KiMc^)) 



■ < ■ 



|det(J^(„))| A'{ui{a))g2{unia)) (1 - d)gi{unia))g2{unia))A{ui{a)) ' 

From Lemma [31 we have that G'{x) > and g'{x) > in M>o- Since G is an 
increasing function, we deduce that 

Gi(u„(a)) 1 



< 



gi{un{a))g2{u„{a)) 2g^(u„(a))' 
Combining the last inequahty with (1341) and (j35p . we obtain 



< 2'l]ada 1 



From Corollary [71 we have that 



(36) 



< 



(l-rf)5^M„(a))A(ui(a)) 
4r/„(r + 1) 



g'2{g~^a)){l-d)a, 
with 0* := 32(5^^(0;*)) (1 — d)a* /2. Hence, defining (5^ in the following way: 

.92(5-'(«))(l-rf)« 



(37) 

we obtain 
(38) 



Sa 



nn < - 



1677„(^* + 1) 



J 



1 

< -. 

- 4 



In particular, for v — </?(«), this bound allows us to consider Ja{u) as a perturbation 
of Jip(a)- More precisely, by a standard perturbation lemma (see, e.g., |OR70[ 
Lemma 2.3.2]) we deduce that Jq(u) is invertible for every u € Bs^{f{a)) and we 
obtain the following upper bound: 



(39) 



4 

< -. 

- 3 



In order to describe our method, we need a sufficient condition for the con- 
vergence of the standard Newton iteration associated to ([7|) for any a £ [a*, a*]. 
Arguing as in |OR701 10.4.2] we deduce the following remark, which in particular 
implies that the Newton iteration under consideration converges. 



Remark 21. Set 5 :— min{(5Q : a G [a*, a*]}. Fix a G [a*, a*] and consider the 
Newton iteration 

, -iF(a,MW) 
< (5, then 
S 



{k>0), 



starting at u^°'> € Ka- If - ^{ce)\\ 

-(^(a)|| 
holds for k > 0. 



< 



Now we can describe our homotopy continuation method. Let ao ■— 
• • • < aN ■= a* be a uniform partition of the interval [a*, a*], with N to be fixed. 
We define an iteration as follows: 

(40) u^^+i) = u^^) - J„,(u('=))-iF(afc,u('=)) (0<fc<A^-l), 

(41) = J,.(u(^+'=))-iF(a*,u(^+^)) (fc>0). 
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In order to see that the iteration (|40l) -(|4T |) yields an approximation of the positive 
solution ^p{a*) of ^ fox A = a* , it is necessary to obtain a condition assuring that 
(PO)) yields an attraction point for the Newton iteration (HT|) . This relies on a 
suitable choice for N , which we now discuss. 

By Theorem [20l we have 

\\^f{ai+i)-'f{ai)\\oo < max{||(^'(a)||oo : a e [a*,a*]} - ttil 

1 AT' 

for < z < — 1, where ki is an upper bound of the condition number independent 
of h. Thus, for N :— \3a*Ki/S~\ + 1 = 0(1), by the previous estimate we obtain 
the following inequality: 

(42) \\(p{a,+i) - ip{a^)\\oo < ^ 

ioi < i < N — 1. Our next result shows that this implies the desired result. 

Lemma 22. Set N := \3a*Ki/5~\+l. Then, for every u'^'^^ with -(/3(a*)||oo < 
5, the point u'^^ defined in is an attraction point for the Newton iteration ^JJ^. 

Proof. By hypothesis, we have Hu'"-* — iy9(a*)||oo < Arguing inductively, suppose 
that II u^*^^ — '■p{ak)\\oo < S holds for a given < A; < A^. By Remark we have that 
M^*^) is an attraction point for the Newton iteration associated to ([7|) for A — Uk- 
Furthermore, Remark [2T] also shows that ||m('^+^^ — </3(afc)||oc < holds. Then 

\\u^^+^^ - ip{ak+i)\\oo < ||w(''+^' -(p(afc)||oo + ||<y5(afc)-<^(afe+i)||oo 

where the inequality ||(/3(afe+i) — '/'(afc)lloo < <5/3 follows by This completes 

the inductive argument and shows in particular that u^^^ is an attraction point for 
the Newton iteration (|4T]). □ 



Next we consider the convergence of (HJ), starting with a point u*^^^ satisfying 
the condition ||u*^^'' — </j(q!*)||oo < a < 5^'. Combining this inequality with ([37]) we 
deduce that u*^^' e Ka* . Furthermore, we see that 

||^t(^+i' -<^(a*)IU = Ik'^' - J.. (u(^')-^F(a*, ^W)-<^(a*)|U 

= (^(^))"^ (J.. (kW)(«W -V(«*)) ~F{a'M''^) + F{a*,v{a'))) 11 

II I I OO 

(43) < ||J..(k(^')-V^(„.)||oc 

|k;i.)(-^c«.(n('^')(^.(^'-^(a*))-F(a*,n('^)) + F(a%^(a*)))|| 
<||J.-(^W)-V^(„.)|U||J-t.)(J..(u('^')-J..(e))||oo||(n('^'-^(a*))||^, 

where ^ is a point in the segment joining the points u^^-' and Lp{a*). Combining 
and ([55)1 we deduce that 

ihi(~+i)-^(a*)iu < |||^;(L.)(^-(«(^))-^-(0)IL'5„. 

< fe<K., 

holds, with c:= {^ria-{6* + 1)) / [g'^^g^^ {a*)){\ — d)a*). By an inductive argument 
we conclude that the iteration (|4T|) is well-defined and converges to the positive 
solution ip{a*) of ([7]) for yl = a*. Furthermore, we conclude that the point u^^^''\ 



24 



E. DRATMAN 



obtained from the point u^^^ above after k steps of the iteration (j4T|) . satisfies the 
estimate 

with c :— 3/4c. Therefore, in order to obtain an e- approximation of ip{a*)^ we have 
to perform log2 log3(3/4ce) steps of the iteration (HH. Summarizing, we have the 
foUowing resuh. 

Lemma 23. Let e > be given. Then, for every u^^'^ £ (R>o)" satisfying the 
condition — (y9(Q:*)||oo < ^, the iteration ^JJ^ is well-defined and the estimate 

_ (^(Q,*)|j_^ <e holds for k >\og2log3{3/ ice). 

Let e > 0. Assmne that we are given it'"^ e (M>o)" such that Hu'-*'-' — </3(a*)||oo < S 
holds. In order to compute an e-approximation of the positive solution (p{a*) of 
dZ]) for yl = a*, we perform N iterations of (|40]) and fco := [log2 log3(3/4c£)] 
iterations of (HTI) . From Lemmas [H] and [53] we conclude that the output u'^+'^o) 
of this procedure satisfies the condition — ip{a*)\\oo < £• Observe that 

the Jacobian matrix Ja{u) is tridiagonal for every a G [a*, a*] and every u S K^. 
Therefore, the solution of a linear system with matrix Ja (u) can be obtained with 
0{n) flops. This implies that each iteration of both (|40| and (|4T|) requires 0{n) 
flops. In conclusion, we have the following result. 

Proposition 24. Let e > and u G (M>o)" with - ip{a^)\\oo <S be given, 
where S is defined as in Remark \21i Then the output of the iteration (^^-(^7p is 
an s- approximation of the positive solution (p{a*) of ^ for A = a* . This iteration 
can be computed with 0(Nn + kt^n) = 0(nlog2 log2(l/£)) flops. 

Finally, we exhibit a starting point u'^^^ G (M>o)" satisfying the condition of 
Proposition [24l Let a* > be a constant independent of h to be determined. We 
study the constant 

5 := minj^Q, : a G [a,, a*]}. 



where 



with 



. ( g',{g-\a))il-d)a , 



C(a 



C2{a) ■.= G-'[G[g-\aC{a)))+ — 
9'2{Ci{a)) 



)a 



2g2{g-\a)/eM)G'{g~\a)le^^y 
C,(<.,:.C-.(c(,-.(^))4), 



and M := g[{Ci{a)). 

Since C(a) > 1, we have that 

7]„ - 2max{5'/(2C2(a)),g^'(2C2(a))a} 

< 2Tnso,{g'l{2C2{a)),g'i{2C2{a))g{G-\G{g-\aC{a)))))}, 

< 2 max{5i' {2G2 (a)) , g'^ {2C2 (a)) <? (C2 (a)) } , 
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As gi and g2 are analytic functions in x = 0, in a neighborhood of € M", we 
obtain the following estimate: 

r]a <2inax{Si{a),S2{a)}{C2ia)Y ^, 

where p is the multiplicity of as a root of gi and Si is an analytic function in 
a; = such that lim^^o Si{a) ^ for i = 1,2. Taking into account that a € (0, a*] 
holds, we conclude that there exists a constant 77* > 0, which depends only on a*, 
with 

for all a G (0, a*]. Moreover, with a similar argument we deduce that there exists 
a constant 1?* > 0, which depends only on a*, such that 



We claim that 

(45) hm = 1+. 

In fact, since we have C{a) > 1 and g~^ is increasing, it follows that 

(46) ^ > 1. 

9 (a) 

On the other hand, there exist ^1 e (^G(g^^(aC{a))),G(g'^{aC{a))) + and 
^2 G (a, aC'(a)) with 

C2{a) ^ ff-i(«g(a))+(G-i)^(6)4 

+ (5-i)'(6)«(5(a) - 1) + (G-i)'(a)4 



< 1 



v{Cia) - 1) a- 



(47) < 



g{g-\a)){Cia)^l) (g(g-^(a)))' 
.9'(.9-H")).9-H«) 2G'(g-i(a))5-H«)' 
Since gi and 172 are analytic functions in x = 0, we see that 

a^o+ g'(g-i(a))g-i(a) 

a^o* 2G'(s-'(a))g-i(a) ' 
Combining these remarks with and (|T7)) we inmediately deduce ps)) . 

Combining (j^J) with (P5|) we conclude that there exists a constant C* > 0, which 
depends only on a* , with 

So.>C*g-\a). 

Therefore, 

(48) (5 = min{(5c, : a e [a,,a*]} > C*g-\l/a,). 
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From Lemma [12] and Lemma [9l we have 
Furthermore, by (|45p . we deduce that 



(49) 



M 



C2{a*)e 



holds for a* > smaU enough. Combining this with (|48l) . we conclude that 
(50) ||u - (p(a,)||oo < C2(a*) - 5~'(a*)/e'*' < 5 

holds for all u G [g-l(Q^*)/e*^ ^2(0*)]". Thus, let a* < a* satisfy dH]). Then, for 
any in the hypercube [(?~-^(a*)/e*^, C2(a*)]", the inequality 

holds. Therefore, applying Proposition 1241 '^^6 obtain our main result. 

Theorem 25. Let e>0 be given. Then we can compute an s- approximation of the 
positive solution of ^ for A ~ a* with 0(nlog2 log2(l/e)) flops. 
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